Imaging diffusional variance by MRI [public] : The role of tensor-valued diffusion encoding and tissue heterogeneity

Diffusion MRI provides a non-invasive probe of tissue microstructure. We recently proposed a novel method for diffusion-weighted imaging, so-called q-space trajectory encoding, that facilitates tensor-valued diffusion encoding. This method grants access to b-tensors with multiple shapes and enables us to probe previously unexplored aspects of the tissue microstructure. Specifically, we can disentangle diffusional heterogeneity that originates from isotropic and anisotropic tissue structures; we call this diffusional variance decomposition (DIVIDE).In Paper I, we investigated the statistical uncertainty of the total diffusional variance in the healthy brain. We found that the statistical power was heterogeneous between brain regions which needs to be taken into account when interpreting results.In Paper II, we showed how spherical tensor encoding can be used to separate the total diffusional variance into its isotropic and anisotropic components. We also performed initial validation of the parameters in phantoms, and demonstrated that the imaging sequence could be implemented on a high-performance clinical MRI system. In Paper III and V, we explored DIVIDE parameters in healthy brain tissue and tumor tissue. In healthy tissue, we found that diffusion anisotropy can be probed on the microscopic scale, and that metrics of anisotropy on the voxel scale are confounded by the orientation coherence of the microscopic structures. In meningioma and glioma tumors, we found a strong association between anisotropic variance and cell eccentricity, and between isotropic variance and variable cell density. In Paper IV, we developed a method to optimize waveforms for tensor-valued diffusion encoding, and in Paper VI we demonstrated that whole-brain DIVIDE is technically feasible at most MRI systems in clinically feasible scan times

Maxwell-compensated design of asymmetric gradient waveforms for tensor-valued diffusion encoding

Purpose: Asymmetric gradient waveforms are attractive for diffusion encoding due to their superior efficiency, however, the asymmetry may cause a residual gradient moment at the end of the encoding. Depending on the experiment setup, this residual moment may cause significant signal bias and image artifacts. The purpose of this study was to develop an asymmetric gradient waveform design for tensor-valued diffusion encoding that is not affected by concomitant gradient. Methods: The Maxwell index was proposed as a scalar invariant that captures the effect of concomitant gradients and was constrained in the numerical optimization to 100 (mT/m)$^2$ms to yield Maxwell-compensated waveforms. The efficacy of this design was tested in an oil phantom, and in a healthy human brain. For reference, waveforms from literature were included in the analysis. Simulations were performed to investigate if the design was valid for a wide range of experiments and if it could predict the signal bias. Results: Maxwell-compensated waveforms showed no signal bias in oil or in the brain. By contrast, several waveforms from literature showed gross signal bias. In the brain, the bias was large enough to markedly affect both signal and parameter maps, and the bias could be accurately predicted by theory. Conclusion: Constraining the Maxwell index in the optimization of asymmetric gradient waveforms yields efficient tensor-valued encoding with concomitant gradients that have a negligible effect on the signal. This waveform design is especially relevant in combination with strong gradients, long encoding times, thick slices, simultaneous multi-slice acquisition and large/oblique FOVs

Diffusional exchange versus microscopic kurtosis from CTI: two conflicting interpretations of the same data

Correlation tensor imaging (CTI) is a new diffusion MRI framework that utilises double diffusion encoding (DDE) to resolve isotropic, anisotropic and microscopic kurtosis sources. Microscopic kurtosis in CTI is provided by the contrast between SDE and parallel DDE signals at the same b-value. Multi-Gaussian exchange (MGE) is a diffusion MRI framework that utilises DDE to measure exchange. The highest exchange sensitivity in MGE is obtained by contrasting SDE and DDE signals at the same b-value. CTI and MGE can thus be applied to analyse the same data but provide conflicting interpretations of that data. We perform Monte Carlo simulations in different geometries with varying levels of exchange to determine which approach is more compatible with the data. Simulations reveal that in all microstructures considered, CTI microscopic kurtosis drastically increases when exchange is introduced. Furthermore, in microstructures that are well-described by the multi-Gaussian assumption, CTI-estimated microscopic kurtosis increases with both the exchange rate and the mixing time, despite fulfilment of the long-mixing-time condition of CTI. Increasing the exchange rate by a factor of 2 positively biases CTI microscopic kurtosis by approximately the same factor. At a modest exchange rate of 10 /s, varying the mixing time from 12 to 100 ms increases CTI microscopic kurtosis by at least a factor of 3. To address this problem, we propose a heuristic approach to combine CTI and MGE to estimate intra-compartmental kurtosis unconfounded by exchange and demonstrate its feasibility using numerical simulations

Gradient waveform design for tensor-valued encoding in diffusion MRI

Diffusion encoding along multiple spatial directions per signal acquisition can be described in terms of a b-tensor. The benefit of tensor-valued diffusion encoding is that it unlocks the "shape of the b-tensor" as a new encoding dimension. By modulating the b-tensor shape, we can control the sensitivity to microscopic diffusion anisotropy which can be used as a contrast mechanism; a feature that is inaccessible by conventional diffusion encoding. Since imaging methods based on tensor-valued diffusion encoding are finding an increasing number of applications we are prompted to highlight the challenge of designing the optimal gradient waveforms for any given application. In this review, we first establish the basic design objectives in creating field gradient waveforms for tensor-valued diffusion MRI. We also survey additional design considerations related to limitations imposed by hardware and physiology, potential confounding effects that cannot be captured by the b-tensor, and artifacts related to the diffusion encoding waveform. Throughout, we discuss the expected compromises and tradeoffs with an aim to establish a more complete understanding of gradient waveform design and its impact on accurate measurements and interpretations of data.Comment: Invited review, submitted in May 2020 to the Journal of Neuroscience Methods. 46 pages, 9 figures, 35 equation

Variability in diffusion kurtosis imaging: Impact on study design, statistical power and interpretation.

Diffusion kurtosis imaging (DKI) is an emerging technique with the potential to quantify properties of tissue microstructure that may not be observable using diffusion tensor imaging (DTI). In order to help design DKI studies and improve interpretation of DKI results, we employed statistical power analysis to characterize three aspects of variability in four DKI parameters; the mean diffusivity, fractional anisotropy, mean kurtosis, and radial kurtosis. First, we quantified the variability in terms of the group size required to obtain a statistical power of 0.9. Second, we investigated the relative contribution of imaging and post-processing noise to the total variance, in order to estimate the benefits of longer scan times versus the inclusion of more subjects. Third, we evaluated the potential benefit of including additional covariates such as the size of the structure when testing for differences in group means. The analysis was performed in three major white matter structures of the brain: the superior cingulum, the corticospinal tract, and the mid-sagittal corpus callosum, extracted using diffusion tensor tractography and DKI data acquired in a healthy cohort. The results showed heterogeneous variability across and within the white matter structures. Thus, the statistical power varies depending on parameter and location, which is important to consider if a pathogenesis pattern is inferred from DKI data. In the data presented, inter-subject differences contributed more than imaging noise to the total variability, making it more efficient to include more subjects rather than extending the scan-time per subject. Finally, strong correlations between DKI parameters and the structure size were found for the cingulum and corpus callosum. Structure size should thus be considered when quantifying DKI parameters, either to control for its potentially confounding effect, or as a means of reducing unexplained variance

Tensor-valued diffusion MRI in under 3 minutes: An initial survey of microscopic anisotropy and tissue heterogeneity in intracranial tumors

Purpose: To evaluate the feasibility of a 3-minute b-tensor encoding protocol for diffusion MRI-based assessment of the microscopic anisotropy and tissue heterogeneity in a wide range of intracranial tumors. Methods: B-tensor encoding was performed in 42 patients with intracranial tumors (gliomas, meningiomas, adenomas, metastases). Microscopic anisotropy and tissue heterogeneity were evaluated by estimating the anisotropic kurtosis ($MK_A$) and isotropic kurtosis ($MK_I$), respectively. An extensive imaging protocol was compared with a faster 3-minute protocol. Results: The fast imaging protocol yielded parameters with characteristics in terms of bias and precision similar to the full protocol. Glioblastomas had lower microscopic anisotropy than meningiomas $(MK_A = 0.29 \pm 0.06$ versus $0.45\pm0.08, p = 0.003)$. Metastases had higher tissue heterogeneity $(MK_I = 0.57\pm0.07)$ than both the glioblastomas $(0.44\pm0.06, p < 0.001)$ and meningiomas $(0.46\pm0.06, p = 0.03)$. Conclusion: Evaluation of the microscopic anisotropy and tissue heterogeneity in intracranial tumor patients is feasible in clinically relevant times frames.Comment: Submitted to Magnetic Resonance in Medicin

Computing and visualising intra-voxel orientation-specific relaxation-diffusion features in the human brain

Diffusion MRI techniques are used widely to study the characteristics of the human brain connectome in vivo. However, to resolve and characterise white matter (WM) fibres in heterogeneous MRI voxels remains a challenging problem typically approached with signal models that rely on prior information and constraints. We have recently introduced a 5D relaxation–diffusion correlation framework wherein multidimensional diffusion encoding strategies are used to acquire data at multiple echo‐times to increase the amount of information encoded into the signal and ease the constraints needed for signal inversion. Nonparametric Monte Carlo inversion of the resulting datasets yields 5D relaxation–diffusion distributions where contributions from different sub‐voxel tissue environments are separated with minimal assumptions on their microscopic properties. Here, we build on the 5D correlation approach to derive fibre‐specific metrics that can be mapped throughout the imaged brain volume. Distribution components ascribed to fibrous tissues are resolved, and subsequently mapped to a dense mesh of overlapping orientation bins to define a smooth orientation distribution function (ODF). Moreover, relaxation and diffusion measures are correlated to each independent ODF coordinate, thereby allowing the estimation of orientation‐specific relaxation rates and diffusivities. The proposed method is tested on a healthy volunteer, where the estimated ODFs were observed to capture major WM tracts, resolve fibre crossings, and, more importantly, inform on the relaxation and diffusion features along with distinct fibre bundles. If combined with fibre‐tracking algorithms, the methodology presented in this work has potential for increasing the depth of characterisation of microstructural properties along individual WM pathways

Impact of velocity- and acceleration-compensated encodings on signal dropout and black-blood state in diffusion-weighted magnetic resonance liver imaging at clinical TEs.

PurposeThe study aims to develop easy-to-implement concomitant field-compensated gradient waveforms with varying velocity-weighting (M1) and acceleration-weighting (M2) levels and to evaluate their efficacy in correcting signal dropouts and preserving the black-blood state in liver diffusion-weighted imaging. Additionally, we seek to determine an optimal degree of compensation that minimizes signal dropouts while maintaining blood signal suppression.MethodsNumerically optimized gradient waveforms were adapted using a novel method that allows for the simultaneous tuning of M1- and M2-weighting by changing only one timing variable. Seven healthy volunteers underwent diffusion-weighted magnetic resonance imaging (DWI) with five diffusion encoding schemes (monopolar, velocity-compensated (M1 = 0), acceleration-compensated (M1 = M2 = 0), 84%-M1-M2-compensated, 67%-M1-M2-compensated) at b-values of 50 and 800 s/mm2 at a constant echo time of 70 ms. Signal dropout correction and apparent diffusion coefficients (ADCs) were quantified using regions of interest in the left and right liver lobe. The blood appearance was evaluated using two five-point Likert scales.ResultsSignal dropout was more pronounced in the left lobe (19%-42% less signal than in the right lobe with monopolar scheme) and best corrected by acceleration-compensation (8%-10% less signal than in the right lobe). The black-blood state was best with monopolar encodings and decreased significantly (p ConclusionAll of the diffusion encodings used in this study demonstrated suitability for routine DWI application. The results indicate that a perfect value for the level of M1-M2-compensation does not exist. However, among the examined encodings, the 84%-M1-M2-compensated encodings provided a suitable tradeoff

Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: Applications in healthy volunteers and in brain tumors

AbstractThe anisotropy of water diffusion in brain tissue is affected by both disease and development. This change can be detected using diffusion MRI and is often quantified by the fractional anisotropy (FA) derived from diffusion tensor imaging (DTI). Although FA is sensitive to anisotropic cell structures, such as axons, it is also sensitive to their orientation dispersion. This is a major limitation to the use of FA as a biomarker for “tissue integrity”, especially in regions of complex microarchitecture. In this work, we seek to circumvent this limitation by disentangling the effects of microscopic diffusion anisotropy from the orientation dispersion.The microscopic fractional anisotropy (μFA) and the order parameter (OP) were calculated from the contrast between signal prepared with directional and isotropic diffusion encoding, where the latter was achieved by magic angle spinning of the q-vector (qMAS). These parameters were quantified in healthy volunteers and in two patients; one patient with meningioma and one with glioblastoma. Finally, we used simulations to elucidate the relation between FA and μFA in various micro-architectures.Generally, μFA was high in the white matter and low in the gray matter. In the white matter, the largest differences between μFA and FA were found in crossing white matter and in interfaces between large white matter tracts, where μFA was high while FA was low. Both tumor types exhibited a low FA, in contrast to the μFA which was high in the meningioma and low in the glioblastoma, indicating that the meningioma contained disordered anisotropic structures, while the glioblastoma did not. This interpretation was confirmed by histological examination.We conclude that FA from DTI reflects both the amount of diffusion anisotropy and orientation dispersion. We suggest that the μFA and OP may complement FA by independently quantifying the microscopic anisotropy and the level of orientation coherence